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Results are presented for the kinetics of domain growth 
of a two-dimensional Ising spin model with competing inter- 
actions quenched from a disordered to a striped phase. The 
domain growth exponent are (3 = 1/2 and /3 = 1/3 for single- 
spin-flip and spin-exchange dynamics, as found in previous 
simulations. However the correlation functions measured in 
the direction parallel and transversal to the stripes are differ- 
ent as suggested by the existence of different interface ener- 
gies between the ground states of the model. In the case of 
single-spin-flip dynamics an anisotropic version of the Ohta- 
Jasnow-Kawasaki theory for the pair scaling function can be 
used to fit our data. 
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The kinetics of phase separation as systems with com- 
peting interactions are quenched below their ordering 
temperature is the subject of much current research in a 
variety of fields in physics and other sciences [ ^ . Some- 
times the equilibrium configuration of these systems is a 
lamellar structure with the ordered phases arranged in 
stripes. Our aim here is to study the growth of striped 
phases in a generalization of the 2D Ising model. In par- 
ticular, we will focus our attention on dynamical proper- 
ties of the correlation functions in the directions parallel 
and transversal to the stripes. 

In the simplest cases the late time ordering process is 
characterized by a single time dependent length R{t) ^ 

representing the average domain size [ ^. This im- 
plies a particular behavior for the correlation function: 
if '^{x, t) is the ordering field the equal time correla- 
tion C[f,t) =< (p{x,t)(f{x + r,!) > has the scahng 
form C{r,t) — f{r/R{t)) where f(z) is a scaling func- 
tion. Our results for the striped phase show the same 
growth exponent in any direction relative to the stripes 
and an explicit dependence of the different longitudinal 
and transversal scaling functions on the microscopic de- 
tails of the system. We will try to understand this dif- 
ference and, in the case of single-spin-flip dynamics, we 
interpret our data by a generalization to anisotropic cases 
of the Ohta-Jasnow-Kawasaki theory [ Q for the scaling 
function. The anisotropic OJK theory is then tested on 
the Ising model with different nearest-neighbor couplings 
for the two square lattice directions. 

Lamellar phases appear in many physical systems. Ex- 
amples are diblock copolymer melts [ |], surfactant- 
oil-water mixtures [ |^, dipolar fluids with long-range 
Coulombic interactions [ Differently than in these sys- 
tems, the striped phase of this paper is not characterized 



by any mesoscopic length - the width of the lamellae is 
always one lattice spacing in the model here considered. 
Neverthless our results may be relevant for systems with 
more realistic lamellar phases. Moreover, from a more 
theoretical point of view, it is interesting to have an ex- 
plicit example of how the scaled correlation functions can 
depend on the details of the system, also in relation with 
recent discussions on this theme [ |^ . 

The model we consider is the well-known Ising version 
of the bidimensional isotropic eight vertex model [ ^, 
with hamiltonian H given by 

- PH ^ Ji ^ SiSj + J2 ^ SiSj + J3 ^ SiSjSkSi, 

<ij> <<ij» [ij.k.l] 

(1) 

where Si are Ising spins on a bidimensional square lattice 
and the sums are respectively on nearest neighbor pairs of 
spins, next to the nearest neighbor pairs and plaquettes. 
Periodic boundary conditions are always assumed. At 
J2 < 0, 1 Ji| < 2| J2I and J3 small a critical line separates 
the paramagnetic phase from a region where four phases 
corresponding to ground states with alternate plus and 
minus columns or rows coexist [ |j. Sudden quenches of 
the system from a completely disordered initial conflg- 
uration to the stripe-ordered phases will be studied by 
numerical simulations both at zero and flnite tempera- 
ture, for different values of the parameters, with heat- 
bath single-spin-flip [ |l^ and spin-exchange dynamics [ 

Before presenting our results it is useful to summa- 
rize the known behavior of system (1) during the phase 
separation process. First consider quenching in the fer- 
romagnetic phase in both d = 2 and d = 3. When all 
the couplings are ferromagnetic the growth exponent is 
/3 = 1/2 or /3 = 1/3 [ |, ]l2| , corresponding respec- 
tively to single-spin-flip and spin-exchange dynamics. 
The situation becomes different for a weak antiferromag- 
netic coupling J2 when, still in the ferromagnetic phase, 
energy barriers oppose to the coarsening of the domains 
and the system does not relax to equilibrium if the tem- 
perature is zero [ In d = 3 the energy barriers are 
proportional to the linear dimension of domains and a 
logarithmic growth is expected [ |l^. Quenching in the 
striped phase in d = 2 have been already studied in [ ^ . 
The simulations of [ ^ show that also in the striped 
phase the average size of domains grows as R{t) ~ 
with (3 — 1/2 or /3 = 1/3 for single-spin-flip and spin- 
exchange dynamics. Motivated by the idea of studying 
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asymptotic correlations in lamellar phases we have ana- 
lyzed again the growth and the dynamical scaling in the 
striped phase of model (1). 




FIG. 1. Typical configuration of model (^) in the scaling 
regime. Grey and white squares represent respectively plus 
and minus spins; black squares represent the interface sites. 
The picture refers to a 100 x 100 square lattice, at zero tem- 
perature, with parameters Ji = 0.1, J2 = —1.0 and J3 = 0, 
after 150 MCS (Monte Carlo Steps per site). The three boxes 
put in evidence different kinds of interfaces. 

A typical configuration of the system during the evolu- 
tion after a quench at T = is given in Fig.l. The config- 
uration is a patchwork of vertical and horizontal striped 
domains. Simulations show that in the case of single- 
spin-flip dynamics there are no energy barriers and the 
system separates also at zero temperature. To monitor 
the domain growth we evaluated the amount of interfaces 
present in the system. To decide if a given site belongs 
or not to an interface, we compare the configuration of 
the system in a neighborhood of the site with four given 
patterns corresponding to the four ground states of the 
model. Then it is possible to define a distance between 
each pattern and the local configuration of the system. 
If this distance is greater than some threshold we say the 
site belongs to an interface, otherwise it is a part of some 
domain [ |l^ . Our results have been shown not to depend 
on the values of the threshold and of the size of the pat- 
terns to be compared. Interfaces identified in this way are 
shown in Fig. 1. In c? = 2 the total length of interfaces 
per unit volume L scales as the inverse of the average 



size of domains. The results of [ |T5[ are confirmed by 
our simulations: we find L ^ t~'^-^ in the case of single- 
spin-flip dynamics and L ^ t^^l'^ in the spin-exchange 
dynamics. The same conclusions have been obtained also 
monitoring the shrinking oi a, N x N square domain of 
one phase immersed in a sea of the three other phases. 
In this case the dipendence on N of the time of shrinking 
can be used to evaluate the growth exponent. 

We now turn to consider the correlation properties of 
growing domains in the scaling regime. We introduce 
the equal time longitudinal and transversal correlation 
functions respectively defined as 

Ce{r,t) = {s{i,j)s{i + e{i,j)r , j + {l-£{i,j))r)) 

Ct{r,t) = ((-l)''s(z,jXz + (l-e(z,j))r-, j+e(i,j)r)) 

(2) 

where e{i,j) is one (zero) if site {i,j) belongs to a hor- 
izontal (vertical) domain and the average is performed 
over all the sites not belonging to interfaces at time t 
and over different stories of the system. These correla- 
tions have the property that for a fixed {i,j) they be- 
come zero outside a given domain. The scaling behavior 
of Cg and Ct for a particular case with single-spin-flip 
dynamics is shown in Fig. 2 where data taken at different 
times have been plotted in terms of the scaling variable 
z = r/^/t. The parameters of the simulations of Fig. 2 
were Ji = 0.1, J2 = —1, J3 = 0. We see that the longitu- 
dinal and the transversal scaling functions fi and ft are 
different, which is a general feature turning out from our 
simulations [ |l^ . Actually the difference between fe and 
ft is very tiny in the case of spin-exchange dynamics, 
but always understandable on the basis of the arguments 
given below. In the following we will show results only 
for simulations with single-spin-flip dynamics. 

To explain the observed difference between fi and ft 
a simple argument can help: a ferromagnetic coupling 
Ji would favour longitudinal with respect to transversal 
correlations while the reverse is true when Ji is negative. 
This is what happens, indeed. The results of simulations 
at Ji = —0.1, J2 = —1, J3 — are the same as those of 
Fig. 2 but with the role of and Ct reversed. Not only, 
fixed J2 and J3, Cg is always greater than Ct for Ji > 
while Ci < Ct for Ji < 0, but the behavior of Ci and Ct 
is very symmetric with respect to the change of the sign 
of Ji. An heuristic argument, based on the existence of 
interfaces with different surface tensions between the 4 
different domains, can explain this symmetry. For sim- 
plicity consider interfaces parallel to the lattice directions 
and J3 = 0. There are three types of these interfaces all 
depicted in the boxes in Fig. 1. On the left of the picture 
a boundary between a vertical and a horizontal domain 
is marked; this interface reduces both longitudinal and 
transversal correlations and it does not matter for our 
argument. The situation is different for the two other 
kind of interfaces marked in Fig.l: the interface in the 
middle reduces only transversal correlations while the one 
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on the right reduces only longitudinal correlations not in- 
fluencing the transversal ones. The T=0 excess energy 
of the interface in the middle is 2 J2 — Ji while the excess 
energy for the interface on the right is 2J2 + Ji. Since 
the excess energy is the driving force for the phase sep- 
aration process, if we assume that the role of these two 
interfaces is the same during this process, we can expect 
a symmetric behavior of the longitudinal and transversal 
correlation functions with respect to the change of sign 
of Ji. 
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FIG. 2. Scaling functions fe{z) and ft{z), z — rjt^'^ in 
the case of single-spin-flip dynamics with T = 0, Ji = 0.1, 
J2 = — 1 and J3 = 0. We averaged over 50 different histories 
on a 400 x 400 lattice. Longitudinal (above) and transverse 
(below) correlations are shown at times 180(«), 220(b), 260(o), 
300(n), 340(A), 380(0), 420(*) and 460(*). The sohd lines 
are the best OJK fits. 

A theoretical prediction for the behavior of the pair 
correlation scaling function in models with non-conserved 
scalar field is given by the OJK theory [ Q . Monte Carlo 
data have been shown to be reproduced by the OJK the- 
ory better than by other theories after an appropriate 
rescaling of the temporal coordinate [ The scaling 
function of the OJK theory is given by 

/(z) = -sin-i[exp(-zV^)] (3) 

TT 

where z = r/t^^^ and D = 8{d — i)/d. This gives the 



Porod linear behavior [ |9[ at small z of the correlation 
function f{z) ^ 1 — az with a = 2^/2/{ny/D). Prac- 
tically, Monte Carlo and theoretical predictions can be 
compared imposing the same Porod behavior. This pro- 
cedure in our simulations gives two different a for the 
longitudinal and the transversal correlations, ae = 0.383 
and at = 0.414 for Ji = 0.1. The above discussed sym- 
metry corresponds to the fact that when Ji = —0.1 the 
best fit with the OJK fmiction is given by — 0.406 
and at — 0.376. 

The above results show that the OJK theory well de- 
scribes the pair scaling function in anisotropic cases if 
a free parameter is used to fit the data in the different 
directions. Since the surface tension is the origin of the 
anisotropic behavior, it is reasonable to check the validity 
of the OJK theory for anisotropic surface tension models 
in the simplest case corresponding to a field model with 
anisotropic kinetic terms. We consider the time depen- 
dent Ginzburg-Landau equation 

where V(ip) is the usual double-well potential. The 
rescaling x ^ x' — \fB^x and y ^ y' = yj Byy would for- 
mally eliminate the anisotropy and would give the usual 
spherically symmetric OJK theory. However the corre- 
lations in terms of the original space coordinates would 
have different Porod laws with ax /ay = By/B^ show- 
ing an anisotropic behavior. Of course, it is not possi- 
ble to simultaneously eliminate by a rescaling both B^ 
and By when higher order derivative terms are present 
in the dynamical equation, so that the scaling function 
in a general anisotropic model is expected to depend in 
some irriducible way on microscopic parameters. 

The anisotropic OJK behavior for the model (4) can 
be tested by studying the dynamical scaling of the Ising 
model with different coupling Jx and Jy in the two lattice 
directions. In Fig. 3 the scaled correlation data in the x 
and the y directions are compared with the OJK scaling 
function. The simulation are at finite temperature with 
Ja; = 2 and Jy — I. Data for the usual Ising model with 
Jx = Jy = ^ are also plotted. The above rescaling ar- 
gument would suggest ax I ay = \J Jy/ Jx- We measure a 
ratio ay/ax — 1.548 which is not far from the expected 
ratio 1.414 and confirms the fact that the OJK theory 
well reproduces the simulation data in anisotropic cases. 
The agreement with the theoretical expectation becomes 
stronger by decreasing the value of the temperature. In 
the case Jx/ Jy — 2 and zero temperature we have mea- 
sured ay/ax — 1.4304. 

To conclude, we have studied domain growth in a spin 
model with four equivalent striped ground states. The 
pair correlation functions are different when measured in 
the direction parallel and transversal to the stripes. An 
explanation of this difference can be given on the basis of 
the different interface energies between the four ground 
states. Our results with single-spin-flip dynamics can 



be analyzed in the context of the Ohta-Jasnow-Kawasaki 
theory which is expected to well reproduce scaling data 
in anisotropic surface tension models. If other quantities 
such as the autocorrelation exponent or other correlation 
functions here not considered depend on the anisotropy 
of the system is a matter for a future study. 
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FIG. 3. Scaling collapse of the correlation functions of the 
Ising model at finite temperature. Data were obtained av- 
eraging over 250 stories in the anisotropic case, Jx = 2 and 
Jy = 1, and over 447 stories in the isotropic case, Jx ~ Jy ~ 1, 
on a 400 x 400 square system. ^From above to below, correla- 
tions along the x-direction in the anisotropic case, correlations 
in the isotropic case, and correlation along the y-direction in 
the anisotropic case are shown at times 350(b), 450 (□) and 
500(A). The solid lines are the best OJK fits. 
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